Thickness dependent Curie temperatures of 
ferromagnetic Heisenberg films 



R. Schiller and W. Nolting 



Humboldt- Universitdt zu Berlin, Institut fur Physik, Invalidenstrafie 110, 

D- 101 15 Berlin, Germany 



Abstract 



We develop a procedure for calculating the magnetic properties of a ferromagnetic 
Heisenberg film with single-ion anisotropy which is valid for arbitrary spin and 
film thickness. Applied to sc(100) and fcc(lOO) films with spin S = | the theory 
yields the layer dependent magnetizations and Curie temperatures of films of various 
thicknesses making it possible to investigate magnetic properties of films at the 
interesting 2D-3D transition. 
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In the past the Heisenberg model in thin films and superlattices has been 
subject to intense theoretical work. Haubenreisser et al. [1] obtained good re- 
sults for the Curie temperatures of thin films [2] introducing an anisotropic 
exchange interaction (2). Shi and Yang [3] calculated the layer-dependent mag- 
netizations of ultra-thin n-layer films with single-ion anisotropy (3) for thick- 
nesses n < 6. Other recent works are aimed at the question of reorientation 
transitions in ferromagnetic films [4] or low-dimensional quantum Heisenberg 
ferromagnets [5]. 

When investigating the temperature dependent magnetic and electronic prop- 
erties of thin local-moment films or at surfaces of real substances it becomes 
desirable to be able to calculate the magnetic properties of the underlying 
Heisenberg model with no restrictions to neither the film thickness n nor the 
spin S of the localized moments. We develop a straigthforward analytical ap- 
proach for the case of Heisenberg film with single-ion anisotropy. 

Considering the Heisenberg model, 
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in a system with film geometry one comes to the conclusion that due to the 
Mermin- Wagner theorem [6] the problem cannot have a solution showing col- 
lective magnetic order at finite temperatures T > 0. 

To steer clear of this obstacle there are two possibilities. First, one can apply a 
decoupling scheme to the Hamiltonian (1) which breaks the Mermin- Wagner 
theorem. The most common example in the case of the Heisenberg model 
would be a mean-field decoupling. For us, the main drawback of the mean-field 
decoupling is its incapability of describing physical properties at the 2D-3D 
transition. 

When choosing a better decoupling approximation to fulfill the Mermin- Wagner 
theorem, the original Heisenberg Hamiltonian (1) has to be extended to break 
the directional symmetry. The most common extensions are the introduction 
of an anisotropic exchange interaction, 

-D^SfS^-D. E (2) 
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and/or the single-ion anisotropy, 
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i iesurf 

In (2) and (3) the first sums run over all lattice sites of the film whereas in 
the second optional terms the summations include positions within the surface 
layers of the film, only, according to a possible variation of the anisotropy in 
the vicinity of the surface. 

Extending the original Heisenberg Hamiltonian (1) by (2) or (3) one can now 
calculate the magnetic properties of films at finite temperatures within a non- 
trivial decoupling scheme. 

For the following we have choosen a single-ion anisotropy which is uniform 
within the whole film leaving us with the total Hamiltonian: 

H = Hf + H A = E J if ( S ia S if3 + S L S jp) + D ° E ( S iaf , ( 4 ) 
ija/3 ia 

where we have considered the case of a film built up by n layers parallel to 
two infinitely extended surfaces. Here, as in the following, greek letters a, 
(3, indicate the layers of the film, while latin letters i, j, number the 
sites within a given layer. Each layer possesses two-dimensional translational 
symmetry. Hence, the thermodynamic average of any site dependent operator 
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A ia depends only on the layer index a: 



(5) 



To derive the layer-dependent magnetizations (S 1 ^) for arbitrary values of the 
spin S of the localized moments we introduce the so-called retarded Callen 
Green function [7]: 



G«t,(B)-«^ ; ^)> E = «^;e° S «S. 
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For the equation of motion of the Callen Green function, 



EG$ a) (E)=h([St,B$ 



+ 



(a) 
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one needs the inhomogenity, 



q+ d(o) 



and the commutators, 



SL H A ] _ = D h (S+S* a + S* a S+) . 



(9) 
(10) 



For the higher Green function on the right hand side of the equation of motion 
(7) resulting from the commutator relationship (9) one can apply the Random 
Phase Approximation (RPA) which has proved to yield reasonable results 
throughout the entire temperature range: 




For the higher Green functions resulting from the commutator (10) this is not 
possible due to the strong on-site correlation of the corresponding operators. 
However, one can look for an acceptable decoupling of the form 

{{S+SL + Sl a S^ B$)) E = $ ia {{Si; B$)) E . (13) 
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As was shown by Lines [8] an appropriate coefficient <&i a = $ a can be found 
for any given function Bj a J = f (Sj^j , which is all we need to know for the 
moment. We will come back to the explicit calculation of the $ Q later. 

Using the relations (8)-(13) and applying a two-dimensional Fourier transform 
introducing the in-plane wavevector k the equation of motion (7) becomes 



(J? (s$ Gf (a) - J? (S' a ) G2f a) ) . (14) 



Writing equation (14) in matrix form one immediately gets the solution by 
simple matrix inversion: 



I (a) 

V o 







71 («0 
In 



(15) 



where I represents the n x n identity matrix and 

h 



(16) 



The local, i.e. layer- dependent, spectral density, = — ^ImCrj^, can then 
be written as a sum of (^-functions and with (15) one gets: 

S£(«0 = E W k )* (E - E,(k)) , (17) 



where E y (k) are the poles of the Green function (15) and XaajO^) are the 
weights of these poles in the diagonal elements of the Green function, G%.f a y 
Both, the poles and the weigths can be calculated e.g. numerically. 

Extending the procedure by Callen [7] from 3D to film structures^ one finds 
an analytical expression for the layer-dependent magnetizations, 

( cz , _ , (1 + Va) 2S+1 (S - <p a ) + <^\S + 1 + <p a ) 

{ a) ' - (1 + ^) 25+1 ' { } 



1 The only pre-condition for the extension is that the spectral density has the 
multipole structure (17) 
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where 



(19) 



Here, iV s is the number atoms in a layer and (3 = The poles and weigths 
in (19) have to be calculated for the special Green function G%? a \ with a = 00. 
In this case the Callen Green function (6) simply becomes: 



n aP _ r a(3 _ // q+, R (0)V 



q+. q- , 



(20) 



and, according to (8), 



V S ) = Va = 2H (S z a ) 



(21) 



Having solved the problem formally we are left with explicitly calculating the 
coefficients $ Q of equation (13). Applying the spectral theorem to (13) for the 
special case of a = one gets, using elementary commutator relations: 



(22) 



We now define the Green function 



D 



/3a 
ji 



S j/3 ;C ia// . 



(23) 



where Ci a is a function of the lattice site. Writing down the equation of motion 
of Dj£* for the limit D -> 0, 



E Dfr(E) = h 



(24) 



and decoupling all the higher Green functions using the RPA one arrives after 
transformation into the two-dimensional k-space at: 



Di a = n 



S 1 , C\ 








\ 



[Sn,C n ]_)j 



{EI 



(25) 



The parameter a had been introduced to derive (19) for arbitrary spin S. 
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where A is a matrix which is independent on the choice of Cj a . Now putting 
C ia in (23) in turn equal to Sf a and to Sf a (2Sf a + ft) and applying the spectral 
theorem to equation (25) one eventually gets the relation: 



rt- q+ 



SjpSf a {2s: a + n) 

Sia> Sf a {2Sf a + ft) 



(26) 



The coefficients & ia are then with (22) given by 



1 i 



Siai Sia(2Sf a + ft) 



rt- o+ 



2((S^)-t^S(S+l) 



(27) 



where, along with commutator relations, the identity 

S&S£ = « 2 S(S + 1) ± " (^) 2 



(28) 



has been used. To avoid the unknown expectation value ((S? a ) 2 ) we apply the 
spectral theorem to the spectral density (17) with a = and get using (21): 



s-s:) = 2H ( S ' a ) j- E E tJ^t ( = 2^ <s;> 

^v s k 7 e i 



(29) 



Hence, with (28) and (29), we get 

(^) 2 ) =ft 2 5(5 + l)- ft (5^(1 + 2^), 



(30) 



and the coefficients can be written in the convienient form 



2h 2 S(S +l)-3h(S z a )(l + 2ip a ) 
(Si) 



(31) 



Together with (31), the equations (15), (16), (18), and (19) represent a closed 
system of equations, which can be solved numerically. 

All the following calculations have been performed for spin S — |, applicable 
to a wide range of interesting rare-earth compounds, and for an exchange in- 
teraction in tight-binding approximation J = O.OleV which is uniform within 
the whole film. The case where the exchange integrals in the vicinity of the 
surfaces are modified has been dealt with by a couple of authors [9]. The 
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Fig. 1. Layer-dependent magnetizations, (S^,), of sc(100) films as a function of 
temperature for various thicknesses n. For all temperatures and film thicknesses the 
{SD increase from the surface layer towards the centre of the films. Inset: Curie 
temperature as a function of thickness of the sc(100) films. 

single-ion ansitropy which plays the mere role of keeping the magnetizations 
at finite temperatures was choosen Dq/J = 0.01. 

Figs. 1 and 2 show the temperature and layer-dependent magnetizations of, 
respectively, simple cubic (sc) and face-centered cubic (fee) films with the sur- 
faces parallel to the (lOO)-planes. For the following Z s means the coordination 
number of the atoms in the surface layers and is the coordination num- 
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Fig. 2. Same as Fig. 1 for fcc(100) films. 

ber in the centre layers of the films. For the case of a monolayer, n — 1, the 
curves for the sc(100) and the fcc(fOO) 'film' are identical, both having the 
same structure. With increasing film thickness the Curie temperatures of the 
films increase. For fcc(100) films the increase in Tq is steeper resulting in the 
limit of thick films in a Curie temperature about twice the value of that of 
the according sc(100) films due to the higher coordination number of the fee 
3D-crystal (Zi,j cc — 12) compared to the sc 3D-crystal (Z^ sc = 6). The larger 
difference between surface and centre layer magnetization of the fcc(100) films 
compared to the sc(100) films can be explained by the lower ratio between Z s 
and Z h (Z s j cc (ioo)/Z b j cc = 8/12 and Z SjSC (ioo)/^6,sc = 5/6). 
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Concluding, we have have shown that the presented approach is a useful and 
straigthforward method for calculating the layer-dependent magnetizations of 
films of various thicknesses and with arbitrary spin S of the localized moments. 
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